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ABSTRACT 

The MMF technique is used to segment the cosmic web as seen in a cosmological N-body simulation 
into wall-like and filament-like structures. We find that the spins and shapes of dark matter haloes are 
significantly correlated with each other and with the orientation of their host structures. The shape 
orientation is such that the halo minor axes tend to lie perpendicular to the host structure, be it a 
wall or filament. The orientation of the halo spin vector is mass dependent. Low mass haloes in walls 
and filaments have a tendency to have their spins oriented within the parent structure, while higher 
mass haloes in filaments have spins that tend to lie perpendicular to the parent structure. 
Subject headings: Cosmology: large-scale structure of Universe - Galaxy formation 



1. INTRODUCTION 

The origin of the angular momentum of galaxies 
and their associated dark matter haloes remains one 
of the most poorly understood subjects in present 
galaxy formation theories despite its relevance in de- 
termining properties such as size and morphological 
type. According to the Tidal Torque Theory (TTT), 
galaxies acquire their angular momentum as a con- 
sequence of the tidal shear produced b y the neigh- 
bouring prim ordial matter distribution • Hovlei 119491 : 
IPeehlesI M969t IDoroshkevkhl 119701 IWhitd I1984T . A 
natural consequence of the TTT is a correlation be 
tween the angular momentu m of haloes and their sur 
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and iLed l|2004D made specific predictions about this. 

Orientation studies based on galaxy cat a logues 
show anti-al ignment llKashikawa fc Okamural 119921 
iNavarro et aT] 120041: iTruiillo et al.ll2006|) . The situation 
in N-body models is less clear: some dark matter N-body 
simulations seem not to detec t any systemic halo spi n 
alignment (iPatiri et alJ l2006t iHev mans et all I2006D. 
while others llLee fc Penll200lMBailin fc Steinmetzll2005l: 
lHatton fc Nininll2001t iFaltenbacher et al]l2002ft present 
evidence that dark matter h aloes are aligned wi th host 
structures. More recently lAltav et alJ l)2006j) found 
a str ong shape alignment of dark haloes in filaments, 
while Bru nino et alJ l)2006j) reported the discovery of a 
systematic alignment effect in an analysis of dark matter 
haloes take n from the very larg e "Millennium" N-body 
simulation (Spring el et aU l20051. We suggest that such 
ambiguities may in large part be a consequence of the 
methods used to identify the larger scale structures that 
host the haloes. 

2. FINDING AND CLASSIFYING STRUCTURE 

Lar ge galaxy surveys such as the 2dF (Coiles set al.l 
I2 QQ1I) or more re cently the Sloan Digital Sky Survey 
l)York et al.ll2000T) have unambiguously revealed an in- 
tricate network of galaxies: the cosmic web. The 
cosmic web can be described as a mixture of three 
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basic morphologies: clusters which arc predominantly 
spherical, elongated associations of gal axies (filaments) 
and large planar structures (walls) ijZeldovichl 119701 
iShandarin fe Zeldovichlll9&^ . 

iBond et al.l l)199fit) emphasized that this weblike pat- 
tern is shaped by the large scale tidal force fields whose 
source is the inhomogeneous ma tter distribution itself 
(see also Ivan de WevgaerU l|200flM . Acccording to TTT 
the same tidal field also generates the angular momen- 
tum of collapsing halos. Thus we would expect the shape 
and angular momentum of cosmic haloes to be correlated 
with one another and with the cosmic web elements in 
which they are embedded. 

Revealing such correlations requires the ability to un- 
ambigiously identify the structural features of the cos- 
mic web. Several methods have in the past been 
used in an attempt to identify and extract the mor- 
phological component s of t h e observed galaxy distri- 
bution iBarrow et all 119851 iBabul fc Starkmanl Il992t 
| Luo fc Vishniad 119951 iStoica et all 120051: iColberg et alJ 
120051 iPimbleti l2005|) with varying degrees of success. 
The results presented in this letter are based on a 
new method, the Multiscal e Morphology Filter: "MMF" 
ijAraeon-Calvo et aD |2006|) . MMF allows us to objec- 
tively segment the cosmic web into its three basic mor- 
phological components by analysing the properties of the 
matter distribution hierarchically. With this morpholog- 
ical characterisation we can isolate specific host environ- 
ments (filaments and walls) for haloes and test predic- 
tions from the TTT in a systematic way. 

The MMF method is a significant advance on other 
similar studies of N-body models. MMF is a technique 
that locates and classifies various structures by exploiting 
localised properties of the inertia tensor of the matter dis- 
tribution on a hierarchy of scales. Since the inertia tensor 
is directly related to the dynamical forces that drive the 
tidal torques, MMF is particularly suited for this inves- 
tigation. The significance of the effects reported here is 
strong, despite the relatively small size of the N-body 
simulation, because of the clear cut MMF environmental 
descriptor. 

3. N-BODY SIMULATIONS AND STRUCTURE 

We ran a cosmological N-body simulation containing 
512 3 equal mass dark matter particles inside a cubic box 
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Fig. 1. — Left panel: Particles inside a sub-box of 37.5 X 75 X 100 h —1 Mpc. For reasons of clarity only a small fraction of the total 
number of particles is shown. Central panel: filaments delineated by a subsample of the particle distribution. At each particle location we 
have plotted the filament vector ep, indicating the direction locally parallel to the filament. Righthand panel: wall particles detected in 
the same sub-box: at each wall particle we plot the wall vector ew Two walls can be clearly delineated: one seen edge-on (dashed outline) 
and one seen face-on (solid outline). 



of 150 h 1 Mpc. using the public vers ion of the parallel 
Tree-PM code Gadget2 l|Snringe]ll20'05T) . We adopted the 
standard cosmological model f2 m = 0.3, f^A = 0.7, h = 
0.7 and <r 8 = 1. The analysis presented here is based 
on the last snapshot at z — 0. The mass per particle is 
2 x 10 9 h _1 M and the softening length was set to 18 
h _1 kpc (comoving) until z = 2 and 6 h _1 kpc (physical) 
after that time. 

3.1. Detecting Filaments and Walls 

The Multiscale Morphology Filter is based on the 
second-order local variations of the density field as 
encoded in the Hessian matrix (d 2 p/dxidxf) for the 
smoothed density field. For a given set of smoothing 
scales we compute the eigenvalues of the Hessian matrix 
at each position on the density field. The density field is 
computed from the particle distribution by means of th e 
DTFE method (|Schaap fe van de Wevgaertjl2000l |2006|) : 
this is natural and self-adaptive, and retains the intricate 
anisotropic and hierarchical features characteristic of the 
cosmic web. DTFE smoothing plays an essential role in 
delineating filaments and walls unambiguously. We use a 
set of morphology filters based on relations between the 
eigenvalues in order to get a measure of local spherical 
symmetry, filamentariness or planarity. The morpholog- 
ical segmentation is performed in order of increasing de- 
grees of freedom in the eigenvalues for each morphology 
(i.e. blobs —> filaments — > walls). The response from the 
morphology filters computed at all scales is integrated 
into a single multiscale response which encodes the mor- 
phological information present in the density field. At 
each stage of the filtering we discard the particles that 
have previously been assigned to a structure and com- 
pute a new density field. 

In a filament, the eigenvectors of the Hessian matrix 
corresponding to the smallest eigenvalue (ep) indicate 
the direction of filament. In walls the largest eigenvalue, 
indicates the perpendicular to the wall (ew). Eigenvec- 
tors are computed from a smoothed version of the density 
field; this avoids small-scale variations in the directions 
assigned to filaments and walls. 

In Figure ^we show a region of the simulation contain- 
ing several filaments and two large walls detected using 



the MMF. The box and its projection were chosen in or- 
der to avoid confusion from projection of many structures 
on top of each other and to show one wall face-on and one 
edge-on. Filaments are clearly delineated like streaming 
lines joining large associations of matter. Walls are more 
difficult to visualise so we also plot their defining eigen- 
vectors. For this projection we can see that the eigenvec- 
tors (ew) defining the wall seen head-on (solid line) are 
pointing towards us, while the ones corresponding to the 
wall seen edge-on are perpendicular to the line-of-sight 
(dashed line). 

3.2. Halo Identification 

Haloes were identified with a somewhat different im- 
plementation of the pub licly available halo finder HOP 
l)Eisenstein fc HuH*1 998). First we identify haloes by run- 
ning hop with the standard parameters and S ou t=&0, 
o~saddie=l20 and 5 peafc =160 for regrouping. Each of 
these haloes is considered a parent candidate which may 
comprise more than one single subhalo. Next for all par- 
ticles we compute densities using a Gaussian window 
with dispersion of 35 h _1 kpc, in order to produce a 
smoothed density field without substructure smaller than 
this kernel. We run hop again but only for particles in- 
side the parent candidates and this time we also provide 
the smoothed densities previously computed as an input 
for hop. The halo identification is performed without 
running regroup. By doing this hop will assign all parti- 
cles to their smoothed local maximum, each group found 
in this stage is a candidate subhalo. For each of these 
new subhaloes we find their center of mass iteratively and 
remove the unbound particles. This "FracHop" method 
allows us to find bound subhaloes inside larger groups 
otherwise identified as single virialised objects. 

Finally, we keep haloes with more than 50 particles and 
less than 5000, a mass range of (1-100) xlO 11 h" 1 M Q . 

3.3. Halo properties 

The angular momentum of a halo containing N parti- 
cles is then defined as: 

N 
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Fig. 2. — A: cartoon illustrating the parameters describing the relative orientations of halos, filaments and walls . B: The probability 
distribution of cos (9: 6 is the angle between the spin vector of a dark matter halos and its host filament. A distinction has been made 
between haloes with masses smaller and larger than 10 12 h — 1 Mq (dashed vs. solid line). The dotted line indicates a uniform distribution 
of halo orientations. The shaded area corresponds to the standard deviation of 1000 random realisations with the same number of galaxies 
as the halo sample. The same prescription is followed to estimate and indicate the standard deviation in the other frames. C: Similar to 
B, for walls: the pdf for the halo spin orientation 8 wrt. the embedding wall. D: The angle cosf between the halo spin direction and the 
host filament as a function of halo mass M. The solid line indicates the average (cos 8) for the halo sample. The dotted line is the expected 
distribution for a randomly oriented sample. G: Similar to D, for walls. Note that the shaded deviation band is wider as a result of the 
lower number of wall haloes. E: The alignment between halo shape and host filament: the average orientation angle (cos 8) of major (solid 
line), medium (dashed line) and minor (dotted line) halo axes relative to the filament. H: Similar to E, for walls. F: The pdf of the angle 
between halo minor axis and halo spin axis for haloes in filaments. The figure distinguishes between low mass (dashed) and high mass 
haloes (solid). I: The pdf of the angle between halo minor axis and halo spin axis, for haloes in walls. Superimposed (shaded area) is the 
same distribution for low mass haloes in filaments taken from F. 



where is the particle's mass, i"; is the distance of each 
particle from the center of the halo, Vi is the peculiar 
velocity of the particle and v the mean velocity of the 
halo with respect to the center of mass. We compute the 
angle between the halo's spin and its assigned Filament 
or Wall, 



are computed by diagonalising the inertia moment tensor 

N 

rriiXiXj, (4) 



= 90 



where 



t>F,w = cos ~ 



J • w 
J || e F , w 



(2) 



(3) 



6f,w being the vector defining the orientation of the fil- 
aments (F) and walls (W). For each halo principal axes 



where the positions of the particles are with respect to 
the center of mass and the sum is over all particles in 
the halo. The orientation of any of the eigenvectors of 
lij relative to the host structure is described by angles 
analogous to the angles 0, 9 W , and Op defined above. 

4. RESULTS 

Our main results are presented in Figure |21 and may 
be summarised as follows: 
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• The orientations of halo spins are significantly cor- 
related with the large scale structure in which they 
are situated. Haloes in walls have spin vectors that 
tend to lie in the plane of the host wall (frame C). 
The strength of the alignment of haloes in walls is 
relatively small (the Lee-Pen correlation parame- 
ter is c = 0.13 ±0.02) but nevertheless statistically 
significant: the K-S probability that the halo ori- 
entations in walls are randomly distributed is less 
than 8.16 x 10~ 5 . We find similar correlations in 
filaments, though of a more subtle character (frame 
B). 

• The orientation of halo spin vectors in filaments 
is mass dependent: low mass halo spins tend to 
lie parallel to their host filament while high mass 
halo spins tend to lie perpendicular to their host 
filament (frames B & D). High mass haloes, M > 
10 12 /i _1 M Q , have a K-S probability less than 5.50 x 
10~ 4 of being randomly orientated, low mass haloes 
less than 1.86 x 10~ 5 (frame B). 

• The principal axes of haloes in filaments are 
strongly correlated with the direction of the fila- 
ments: the minor axis tends to be directed perpen- 
dicular to the filament, the effect being strongest 
for the larger masses (frame E) 

• The principal axes of haloes in walls are strongly 
correlated with the perpendicular to walls: the mi- 
nor axis tends to be lie perpendicular to the wall 
while the other axes tend to lie in the local plane 
of the wall (frame H). 

• While both the halo minor axis and the spin are 
correlated with the megaparsec scale structure, 
the distribution of the angle between these two is 
skewed and biased towards the surrounding large 
scale structure. This slight skewness may bias the 
estimate of the mean angle the spin makes with the 
host structure (frames F and I). 



5. CONCLUDING REMARKS 

The significance of the alignments we have found, given 
the relatively small sample of haloes and small simulation 
volume, emphasises the benefits of having a good def- 
inition of environmental structures when studying halo 
properties and their relation with the cosmic web. 

The results r eported here are in ac cord with results re- 
ported by e.g. IBrunino et al.l l|2006j) . We find that halo 
spin vectors lie within the host structure and halo minor 
axes point out of the host structure as previousl y sug- 
gested by ijPatiri et alJl2006t IBrunino et all 12006). The 
major advance of our study is the objective and multi- 
scale character of the structure identification by means 
of the MMF. Earlier studies at best resor ted to heuristic 
means of delineatin g filaments or walls l|PimbletJ 120051 
iColberg et all 120051: Brunino et alJl2006|) . Some earlier 
studies have measured the distance of a halo from a local 
minimum in the large scale density distribution. On the 
assumption that the surfaces of such "spherical" voids 
are to be identified with the large scale structures hosting 
the haloes one may attempt to find signatures of spin and 
shape alignments. Voids are not generally spherical, nor 
is it always clear which of several voids a halo is related 
to. Moreover, using voids as defining the environment of 
a halo does not in itself allow for the important distinc- 
tion between walls and filaments. Alignment detections 
should therefore be seen as residuals of genuine physical 
alignments with walls or filaments. 

This letter poses a number questions for future investi- 
gation. The mass segregation of alignments in filaments 
is yet to be understood. Detailed merging trees may give 
clues to the possible non-linear processes responsible for 
this effect. It will of course be interesting to use MMF to 
define environments in catalogues of galaxies such as 2dF 
and SSDS: this sh ould provide an impo rtant refinement 
of earlier results of lTruiillo et alJ l|2006j) and since it will 
unambiguously localise galaxies as either being in walls 
or in filaments. 

We thank Pablo Araya for the N-body simulations used 
in this work and gratefully acknowledge discussions with 
Ignacio Trujillo during the "Bernard60" conference in Va- 
lencia. 
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